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When a medium composed of microscopic elements is subjected to a iiigii intensity field, the 
individual behaviors of microscopic elements can become chaotic. In such cases it is important to 
consider the effects of this irregularity at microscopical level onto the macroscopic behavior of the 
medium. We show that the macroscopic field produced by a large group of chaotic scatterers can 
remain regular, due to the partial or complete phase coherence of the scattering elements and the 
incoherence of the chaotic components of their responses. Thus when only macroscopic fields are 
observed, one may be unaware of chaotic microscopical motion, as it appears to be hidden from 
the observer. The coupling among the elements may lead to partial chaos synchronization, which 
exposes the chaotic nature of the system making the oscillations of macroscopic fields more irregular. 



The problem of wave scattering and dispersion is one 
of the basic problems in many branches of physics. In 
most physical problems scattering and dispersion are ex- 
plained by the interaction of the wave field with particles 
that respond linearly or weakly nonlinearly to the exter- 
nal forcing. However over the last few years problems 
appeared where external fields are so powerful that they 
drive the individual scatterers into the chaotic regime. 
Examples of such problems are the chaotic excitation of 
atoms and molecules by powerful electro-magnetic fields 
and chaotic pulsations of cavitating bubbles in power- 
ful ultrasonic fields |^ . In such cases the problem of scat- 
tering and dispersion becomes more complicated. Even 
if the behavior of a microscopic element is chaotic, it is 
not obvious what manifestation this would have at the 
macroscopic level. The study of this issue is important 
from both theoretical and experimental prospectives. 

Scattering elements can typically be considered as pas- 
sive damped oscillators. When the amplitude of the scat- 
tered field is small compared to the external field, the 
effect of the scattered field upon an individual element is 
negligible and the whole system becomes an ensemble of 
uncoupled damped oscillators driven by the same exter- 
nal force. Beyond some amplitude level the nonlinearity 
of the scatterers becomes essential and can lead to chao- 
tization of their oscillations. In many cases chaotic oscil- 
lations occurring in a single damped oscillator under the 
action of a periodic external force turn out to be in some 
sense phase locked to the driving force As the result, 
chaotic oscillations in such systems contain components 
that are coherent with the external driving. If the macro- 
scopic field is produced by a large group of such micro- 
scopic elements, these components add coherently, while 
the chaotic components of signals add non-coherently. 
The resulting macroscopic field becomes nearly periodic. 

When the coupling through the scattered field is taken 
into account, an unusual situation occurs. Sufficiently 
strong coupling can lead to partial or complete synchro- 
nization of chaotic oscillations in individual elements [Q . 
In such case chaotic components of these oscillations be- 
gin to add coherently and the macroscopic field becomes 
more chaotic. Thus the synchronization phenomenon. 



which is normally associated with the onset of a more 
regular behavior in the system, has the opposite effect 
on the observed macroscopic quantities, making their os- 
cillations less regular. 

To better understand this phenomenon, let us consider 
as an example the following system: 
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These equations describe damped motion of N par- 
ticles in the potential U{x) = a/(2a;^) — 1/x and sub- 
ject to the periodic external force. The last term is to 
account for possible synchronization effects due to the 
mutual coupling through the macroscopic field. Here we 
assume for simplicity that the scattering region is much 
smaller than the wavelength of the external field and ne- 
glect the delays in the coupling term. The values of the 
coupling coefficients depend on the geometry of the prob- 
lem. We consider the case where the scattered field is 
proportional to the derivative of the state of a scattering 
element, yj{t) = dxj{t)/dt. In the following discussion 
we shall keep i' = 0.4, a = 0.75, A = 0.45, and = 1.3. 
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FIG. 1. The chaotic attractor in one oscillator 
coupling. 

With these values of parameters and with Ki j = for 
all i and j, system (|l|) oscillates in a chaotic regime, as 
illustrated in Figj^. 

The response of the ensemble, Y{t), observed far from 
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the scattering region is proportional to the sum of fields 
from individual scatterers: 



Y{t) 



N 



Assuming stationarity of yj{t), the mean temporal 
power spectral density P{u!) of the scattered field can 
be calculated using j|] 



Pico) 



C{t) coslot dr, 



where C(r) is the auto-correlation function for Y{t): 



N 



C{t) =< Y{t)Y(t + T) >= 



The last sum contains auto-correlation terms (j = k) 
and cross-correlation terms (j 7^ k). 

Let us first consider the case when Kij = for all i and 
j. Then all oscillators (|l|) evolve on the same attractor so 
that < yj{t)yj{t + T) Cq{t), and < yj{t)yk{t + T) 
Cx{t) for i =^ j, independent of j and k. Thus the 
autocorrelation of the macroscopic field can be written 
as 

C{t) = N {Coir) - Cx{t)) + N^Cx{t), 

and the power spectral density can be written as 

PiLo)^N{Pi{uj)-Px{Lo)) + N^PxH, (2) 

where Pi{uj) is the power spectral density of a single re- 
sponse, and Px{lo) is the cross spectral density. 

The autocorrelation, Co(t), the cross-correlation 
Cx{t) and the difference between the two, Cc(t) — 
Cq{t) — Cx{t) are shown in Fig.^. We see that Cc(t) 
decays at large |t|, which means that -Pi(w) — Px(w) is 
the continuous component of the power spectrum. More 
interestingly, we observe that Cx (t) is a purely periodic 
function, meaning Px(^) is the discrete component of 
the spectrum. 

Equation || indicates that as the number of scatter- 
ing elements increases, the coherent component of the 
spectrum increases quadratically with the number of el- 
ements, while the continuous spectrum component in- 
creases linearly, as it typically happens with non-coherent 
signals. Therefore, with a large number of chaotic micro- 
scopic elements, the macroscopic response of the system 
is highly regular, nearly periodic. Thus the chaotic na- 
ture of the system is hidden from the macroscopic ob- 
server. This is illustrated in Fig.|| and Fig.^. The faster 
growth of the discrete spectral component, compared to 
the continuous chaotic component, is evident in Fig.^ A 
closer look at a proper projection of this figure reveals 
good agreement with the prediction of (ph: the slope for 



the discrete part of the spectrum is twice larger than for 
the continuous part. 
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FIG. 2. Autocorrelation of the chaotic response of oscilla- 
tor (^) without coupling (top); cross-correlation of responses 
of two oscillators (middle); and the difference between these 
correlations (bottom). 

In order to explain the mentioned features of the cor- 
relation functions, let us study the phase properties of 
individual responses. System (|l]) with zero coupling be- 
longs to the class of two-dimensional systems with pe- 
riodic forcing, dx/dt — F(x, i), where x G and 
F(x, i) — F{x,t + T), with strictly negative divergence 
of the vector field x: divF(x, t) < for any x. By Liou- 
ville's formula it can be shown that for such systems the 
Poincare map x(i) — > x(t + T) is area-contracting and 
therefore cannot contain invariant circles Q. For the 
non-autonomous system x(t) this means that it cannot 
posses quasiperiodic motions corresponding to trajecto- 
ries on invariant tori. All motions in this system are 
either periodic motions phase locked with the external 
forcing or chaotic motions. If present, chaotic trajecto- 
ries will wonder in the vicinity of the skeleton composed 
of periodic orbits, which arc phase locked to the driving 
force 01 . The effect of this phase locking can be seen if 
we introduce the phase of chaotic oscillations as the an- 
gle variable in the angle-amplitude representation of the 
Hilbert transform, hj{t) 



h,{t)^y,{t) + -FN. 



Vj{t) 



t 



dr. 



The phase mismatch between the response and the 
driving force, ^j(t) — (t>j(t) — fit, is shown in Fig||. 
We observe that max(<i>j(t)) — minj(<i>(t)) < tt. Thus 
the chaotic response oscillations are loosely phase locked 
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with the periodic driving and the fluctuations of the rel- 
ative phase are less than half-a-period. Due to this, the 
cross-correlation of two oscillators driven by the same ex- 
ternal field is purely periodic because the effect of chaos 
averages to zero, while the component coherent with the 
periodic driving does not. When t is large, due to the 
divergence of trajectories in the phase space of chaotic 
systems and the loss of information about the initial con- 
dition, the difference between computing autocorrelation 
for a single oscillator and computing cross-correlation for 
signals from two different oscillators disappears. This 
explains why Co(t) becomes purely periodic, and Cc(t) 
decays to zero, see Fig.|[ 
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FIG. 3. The mean field response of 4096 oscillators. The 
response of a single element is shown with a dotted line for 
comparison. 




FIG. 4. The mean power spectral density of the macro- 
scopic response as a function of the number of oscillators (in 
logarithmic scale). 

Let us now consider what happens when the coupling 
among the elements is taken into account. For simplic- 
ity we assume that Kij = K. Then the system (|l|) has 
the solution where Xi{t) — X{t) for all i, which corre- 
sponds to identically synchronized chaotic oscillations in 
the ensemble. If the coupling among the elements is suffi- 



ciently strong to stabilize this solution, clearly, the mean 
field response is proportional to that of a single element 
and is chaotic. Thus, although as synchronization sets in 
the dynamics of the system becomes less chaotic, with a 
lower dimension of the attractor, the mean field response 
becomes more irregular. 
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FIG. 5. The detrended phase of the response of an individ- 
ual oscillator. 
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FIG. 6. The ratio of the power in the continuous compo- 
nent of the spectrum to the total power. The solid line cor- 
responds to increasing the coupling coefficient from zero, the 
dashed line, to decreasing it. 

Not quite so obvious is the effect of weak synchronizing 
coupling. When the coupling is weak, the state of com- 
plete synchronization may not be achieved. To make the 
matter more complicated, the variations of coupling lead 
to changing dynamics of the entire systems. In addition 
to chaotic attractors for some values of coupling periodic 
or quasiperiodic stated can become stable. Nevertheless, 
the de-regularizing effects can be seen even for very small 
values of the coupling. This is illustrated in Fig.^ which 
shows the ratio of the power in the continuous compo- 
nent of the spectrum to the total power as a function of 
coupling. We see that for weak coupling we observe a 
sharp transition from a very regular state where most of 
the power is in the periodic components to a less regular 
state with a significant power in the continuous part of 
the spectrum. A hysteresis is observed near the transition 
point. As the coupling increases further, the dynamics of 
the system changes and the mean field response becomes 
more regular again. At larger values of coupling, 0.2, 
the system begins to approach completely synchronized 
state. Kt K = 0.2 the cross-correlation of of responses 
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of two oscillators looks almost identical to the autocorre- 
lation of a single oscillator, and the oscillations become 
quite irregular. 

In conclusion, we showed that when the only observed 
quantities are obtained by averaging the responses from 
many chaotic oscillators to external fields, these quan- 
tities can remain periodic. This presents a difficulty in 
observing chaos in such systems. Fig.^ illustrates the 
effect of averaging onto the bifurcation diagram of our 
example system. We see that the period doubling cas- 
cade and the chaotic regime, evident in the bifurcation 
diagram of a single element, are hardly visible in the bi- 
furcation diagram of the mean field. Under experimental 
conditions the bifurcation sequence and the transition to 
chaos can easily be obscured by measurement noise [||. 
The macroscopic chaotic oscillations may arise in such 
systems as a result of mutual coupling between the ele- 
ments. When this coupling is strong chaotic oscillations 
in individual elements may become synchronized, which 
causes non-periodic components of individual responses 
add coherently. As the coupling strength varied, the tran- 
sition between regimes characterized by different degrees 
of regularity of the mean field response can occur in a 
non-trivial way, for example, exhibit hysteresis. 
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FIG. 7. Bifurcation diagrams for a single oscillator and for 
the mean field of 4086 elements (bold line) . The diagram was 
created using the Poincare map with the period of the driving. 
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Partial coherence and effects similar to those discussed 
in this communication can also occur in ensembles of 
mean field coupled chaotic generators. In such systems 
the external field may not be necessary to achieve phase 
coherence, which can arise spontaneously due to a global 
mean field coupling in the system [^ . 
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